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ABSTRACT 

We have compared Monte Carlo photoionization models of HII regions with a 
uniform density distribution with models with the same central stars and chemical 
compositions but with 3-D hierarchical density distribution consisting of clumps within 
clumps, on a four-tier scheme. The purpose is to compare the abundances of He, N, O, 
Ne, and S obtained by standard analyses (emission line strengths and measured mean 
temperatures from [OHI] and [Nil]) to the abundances in our models. We consider 
stellar temperatures in the range 37.5 - 45kK and ionizing photon luminosities from 

1048 _ IQSl g-i. 

Clumped models have different ionic abundances than uniform. For hot stars, 
(He°/H+) is 2 - 3%, much larger than with uniform models. This amount of He° is 
independent of metallicity and so impacts the determination of the primordial abun¬ 
dance of He. The total abundances of O, Ne, and S obtained by the usual methods of 
analysis, using r([OHI]) for high stages of ionization and r([NH]) for low, are about 
as accurate for clumped models as for uniform and within ~20% of the true values. 
If T([OHI]) is used for analyzing all ions, the derived (0/H) is ~ 40 - 60% too large 
for cool stars but is good for hot stars. Uniform models have similar errors, so the 
clumping does not change the accuracy of abundance analysis. 

The physical causes of the ionic abundance errors are present in real nebulae. 
In clumped models, helium ionizing radiation from zones of high ionization (low He° 
and low UV opacity) can penetrate nearby regions near the edge of the ionized zone. 
This effect allows He° to absorb more stellar photons than in uniform or radially 
symmetrical geometries. In turn, these absorptions compete with O'*', etc., for those 
energetic stellar photons. 

Key words: radiative transfer — ISM — H H regions — cosmology 


1 INTRODUCTION 

The interstellar medium (ISM) has filamentary structure 
over scales ranging from hundreds of parsecs through AUs. 
This paper will analyze some of the properties of H II re¬ 
gions formed by ionization of the clumpy interstellar gas 
and compare them to those produced by the same code for 
models with a spatia lly uniform densit y. The Monte Carlo 
photoionization code iWood et al .l2004^ treats gas with den¬ 
sities that have a power-law or scale-independent character 
over about an order of magnitude in spatial sizes (see §2). 
We refer to these hierarchical models as “clumped”. 

Our goal is to isolate the effects of clumping by com¬ 
paring clumped models with smooth density distributions, 
while keeping other important parameters (e.g., the stel¬ 


lar atmospheres and luminosities that strongly influence the 
nebular models, and the nebular composition) constant. We 
choose the other parameters rather arbitrarily and do not 
try to compare our models with real objects. Several types 
of questions motivate us: 

(a.) For a fixed exciting star and mean nebular density, 
how are the ionization and temperature structures of the 
resulting ionized nebulae affected? What are the physical 
causes of the differences? How do the differences vary as the 
nebular parameters (e.g., exciting star properties, geometry 
of the clumps, and the mean gas density) are changed? 

(b.) How accurately are model ionic abundances deter¬ 
mined from the emission lines, using techniques similar to 
those commonly used to analyze observations of real nebu- 
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lae? What errors occur in the ionic abundances derived from 
ratios of fluxes of emission lines, using mean nebular electron 
temperatures, Te, estimated from diagnostic line ratios? 

We assume idealized conditions, neglecting observa¬ 
tional errors and reddening of emission lines by interstellar 
dust either within or exterior to the HII region. Thus, our 
abundance errors are lower limits to those expected for real 
nebulae. We consider the effects of dust in the transfer of 
Lyman continuum radiation within the ionized gas. 

Our assumptions and the limits of parameter space that 
we explored are given in §2. In §3 we consider the physical 
differences among the clumped models and their differences 
from uniform models. In §4, we discuss the relevance of our 
models to real nebulae and the uncertainties of abundances 
within them. 

We do not address the important issue of differences in 
abundances of many ions derived from recombination lines 
as opposed to collisionally excited lines. A common expla¬ 
nation is that there are “anomalous” te mperature variations 
withi n the ionized plasma (see, e.g., lGarcia-R,oias et alJ 
I 2 OO 4 I . and references therein) , although the ph enomenon 
seems more complicated (e.g., iTsamis et al.l20f)^ . We have 
used only standard physics, and our abundances from recom¬ 
bination lines would be almost the same as from collisionally 
excited lines. 


2 ASSUMPTIONS USED IN THE MODELS 

There are many free parameters that affect the spectrum 
and abundances derived from model H II regions, and we 
consider only a subset of possible parameter space. We make 
the following assumptions: 

1. The models are in a steady state, thereby neglecting 
the dynamics occurring near the outer H'^ —> H ° ionization 
front , where probably a D-type front (see, e.g.. lOsterbrockI 
Il989^ is propagating into the surrounding neutral material. 
This is the standard assumption in H II region models. How¬ 
ever, the conditions near the ionization front are important 
for N° and 0°, which are sensitive to the temperature and 
ionization structure at the outer edge of the ionized region. 

2. We chose the composition used for workshops 

for testing model H _ II regions iFerland et al.l I199RI 

IPeauignot. D.. et ahlliioih : H, He, C, N, O, Ne, and S 
= 1, 0.1, 2.2 4.010“®, 3.3 10"'^, 9.0 10"®, and 5.010“®, 

respectively, by number. These abundances are rather arbi¬ 
trary, but are adequate for investigating the important phys¬ 
ical effects introduced by clumping. We will refer to elemen¬ 
tal abundances as chemical symbols without superscripts, so 
that (0/H) is the abundance ratio of O, including all stages 
of ionization, relative to H. 

3. Dust in HII regions and within the diffuse ionized 
gas (DIG) needs careful consideration. We will show below 
that carbonaceous grain material has probably largely been 
converted to the gas phase. We will find that the effects of 
carbon-depleted dust are not very significant, but that with 
other assumptions dust could become dominant in deter¬ 
mining the structure of the ionized region. 

Dust affects the emitted nebular spectrum by reddening 
the nebular radiation and also by absorbing ionizing radia¬ 
tion. Since corrections for reddening are routinely made, we 
assume that line strengths are determined perfectly. This as¬ 


sumption is very optimistic when we consider the spatially 
unresolved emission line fluxes rather than the pixel-by-pixel 
line intensities (see §3.2). 

If there were no destruction of dust within the ionized 
zone, the absorption of ionizing radiation by dust would be 
hugely important, while dust scattering, strongly directed 
forward, has almost no effect. The absorption cross section 
of graphitic carbon, closely related to absorption by poly- 
cylic aromatic hydrocarbons (PAHs), peaks at about 17 eV 
with a cross se ction of 3 x 10“^^ cm^ (H ato m)“^ according 
to the model of lWeingartner fc Draing ll200l|). wit h the cross 
sections conveniently tabulated iri T^ainel ^2004^ . This ex¬ 
ceeds the opacity of H if (H®/H)^; 0.001, which occurs over 
most of the volume of the nebula. Dust with this carbona¬ 
ceous cross sections absorbs ~90% of the ionizing radiation 
over our entire range of models, for all geometries. By con¬ 
trast, the silicate component of dust is relatively gray in the 
ionizing UV, with an opacity of roughly 1.3 x 10“^^ cm^ 
H“^, only 5% of the peak of carbonaceous dust. 

The abundance of dust in H H regions can be estimated 
by means of the gas-phase abundance s of refractory ele¬ 
ments (e.g., C, Al, Mg, Ca, or Fe; see ISavaee fc Sembadil 
I 199ftl that are mostly locked up in dust in the diffuse 
ISM. Abundances found in bright HII regions are sig¬ 
nificantly different, especially for C. In the diffuse neu¬ 
tral ISM, the CII abundances derived from absorption to¬ 
wards 12 B stars is C/H = 1 60T17 ppM and an upper 
limit of 108 ppM for another llSofla et alJl2003) . In vari- 
ous HII regions, gas phase C/H varies from 330 to 590 ppM 
j|Esl^ari_gtaljj2004a|jEstel^^ alJ fT999aL lEsteban et alJ 
f^99bl ~lGarci^Roias~e^ni200^. Solar (C/H ) = 245T30 
ppM if^phmr^^in^OOSnHolwegeil |2^^. Interstellar 
abundances have probably been significantly enhanced since 
the fo rmation of the Sun, p e rhaps dex(0.2) to ~400 ppM 
(e.g., lAckerm an etjajJjgOT^^JPage^ ^ Tautvaisiend ^9^. 
The Orion Bar dSellCTeiretSTIl^Onl) . an ionization front 
seen edge-on, shows 3.3 /rm aromatic hydrocarbon emission 
arising from a rather thin zone between the Pq and Bra 
recombination lines of H and emission lines of H 2 . Thus, 
several observations require that bright HII regions have 
most or almost all of their carbon in the gas phase. Usually 
HII regions have strong near-infrared emission bands, but 
it may be that the carriers are actually in the neutral gas 
adjacent to the ionized. 

In some HII regions, silicaceous dust is present within 
the ionized gas. In M17, only about 15% of the Fe is 
in the gas phas e llEste ban et al1ll999bh . Within the DIG, 
iHowk fc Savag^ i 199fll) found that AlIII/SlII is much less 
than solar. Both ions are found only in the gas in which H 
is ionized, so Al is depleted onto grains in the DIG. Further 
direct evidence for dust within the DIG is given by the cor¬ 
relation of the diffuse H a with 100 /rm emissi on and the lack 
of correlation with HI iLagache et a il l 200(Tl) . 

In summary, it seems that the carbonaceous component 
of dust has largely been vaporized in HII regions while the 
silicaceous grains may largely remain intact. Perhaps pro¬ 
cesses that convert carbon in grains to the gas phase act 
very rapidly in ionized gas, so that the C in dust (including 
PAHs and the “very small grains,” carbonaceous in charac¬ 
ter) does not absorb ionizing radiation when H is mostly H^. 
Pure silicate dust has cross sections so low and so weakly 
dependent on wavelength that its presence is not crucial to 
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the results of this paper. We ran models that confirmed this 
statements, but for most models in this paper we neglected 
dust and also heating by PAHs. 

4. Our algorithm for clumping was not derived from any 

model of turbulence, but rather from hierarchically clumped 
m odels of the density distribut i ons in the manne r described 
in iMathis. Whitney, fc WoodI following lElmegre^ 

<199711 . Points were placed within a cube as described be¬ 
low. The cube was partitioned into 65® cells, and the final gas 
density within each of the cells was proportional to the num¬ 
ber of points falling within the cell. The exciting star was 
placed at the center of the cube. Initially 32 points were cast 
randomly throughout the cube. Surrounding each of these 
another 32 points were cast randomly within a distance de- 
termined by the “fractal dimens ion”. D, of the model (see 
iMathis. wiiitnev. fc WoodI I 2 OO 2 I for details). Around each 
of these points another 32 were nested, and finally 32 more 
around each of the preceding round of casting. Thus, the 
structure, if it were continued indefinitely, would be con¬ 
sidered “fractal”, but the power spectrum is a power law 
only over about a decade in sizes of projected densities. The 
mean density was taken to be n.(H) = 100 cm“®, with a 
fraction /smooth in a uniform distribution and the rest in the 
clumped cells. The /smooth represents an average over a still 
finer distribution of subclumps, not spatially resolved in our 
simulation, that extend down to very small scales. 

For this paper we chose /smooth = 0.15, so the mini¬ 
mum density of a cell is 15 cm“®. The maximum density 
within these models was 2700 cm“®, with a spread of ±15% 
arising from different initial seeds. For comparison to the 
clumped geometries, uniform models were computed with 
the same stellar luminosity and n(H) = 100 cm“®. The uni¬ 
form models fitted more compactly into the grid cube and 
consequently could employ a mesh 2.7 times finer than the 
clumped models. The abundance patterns found with the 
coarse and finer resolution uniform models are in close agree¬ 
ment except for 0° (the coarse being 1.5 - 2.2 times larger). 

In order to test the effects of allowing radiation to es¬ 
cape from our clumped models through vacuum instead of 
unresolved cells, we ran models with /smooth = 0. About 10% 
of the ionizing photons escaped the simulation, as opposed 
to ~1.5% for the /smooth = 0.15 models. 

The distribution of the density in the clumped models is 
determined by three parameters: (a) /smooth, (6) the initial 
seed for the random number generator, and (c) the “fractal 
dimension”, D, that determines the power spectrum of col¬ 
umn densities as projected from an assumed cube onto the 
plane of the sky. We considered D — 2.6 and 2.9, and av¬ 
eraged all clumped models among the same five arbitrarily 
chosen initial seeds. Thus, we have an indication of the dis¬ 
persion of various physical quantities among clumpy models 
with the same fractal dimension. 

5. The H-ionizing photon luminosity of the exciting 
star in units of 10^® photons s“^, L49, was increased from 
0.1 to 100, in steps of factors of ten. For uniform models 
this is equivalent to changing the “ionization parameter”, 
U = 10"'® L49/(47ri7| riec), by factors of lO'^® (Rs is the 
Stromgren radius and the electron density.) A large U 
implies large abundances of high stages of ionization and 
a relatively sharp increase of (H®/H) near the edge of the 
Stromgren sphere. In clumped models the ionized region is 
not spherical and there is no constant value of rie, but L49 


is an equivalent parameter to explore the effects represented 
by U in uniform models. For instance, our results are appli¬ 
cable to superstar clusters: ±49 ~ 10 ® — 10 ^, rie ~ 10 ® — 10 ^ 
cm“® in starburst galaxies. 

6 . The dielectronic recombination coefficients for sul¬ 

fur are poorly known. Unfortunately , they can be imp ortant 
for the ionization of S. We followed lAli et akl using 

averages of various ions of C, N, and O. We tested the impor¬ 
tance of these assumptions by calculating some models with 
all dielectronic recombinations for S being zero. Of course 
this is an extreme assumption, since a lack of knowledge of 
the values of coefficients does not imply that they are zero. 
We discuss the results in H3.5.4I 

The cube for clumped models was large enough so that 
< 1 % of the ionizing radiation escaped through the edges. 
The outer boundaries of the ionized regions were very irreg¬ 
ular in shape, and for all stellar temperatures there are a 
few regions of neutral H and He embedded within the ion¬ 
ized gas, with the gas that they shadow from the star ionized 
by recombination radiation. 

7. Our cell sizes are large in comparison to the mean free 
paths of ionizing photons when H is almost neutral within 
the cell. Our code derives a mean ionization for the cell 
based on the radiation entering it without considering any 
substructure of the ionization within the cell. Our analy¬ 
sis attributes that ionization to the entire mass of the cell, 
while only a fraction of the cell’s mass would actually be in¬ 
volved. Near the edge, a few photons can produce ionizations 
that the code assumes are spread uniformly over the entire 
cell, with a correspondingly low electron density. Therefore, 
recombinations are underestimated, and too much mass of 
ionized material is derived. The principal ions for which this 
effect occurs are H’*’, N®, and O®. In our analysis we con¬ 
sidered only cells in which the code derived (H'^/H) less 
than H^in, a parameter taken to be 0.25 for the results in 
this paper. This choice provides roughly an optical depth of 
unity across our cells, for a typical photon that can reach 
the edge. We have investigated how these results depend the 
Hmin- Uniform models are hardly affected because the area 
of their almost-neutral H zone is much smaller than clumped 
models. These show a slow decline in (He''"/H''") as is 
increased from 0.1 to 0.4, and a more rapid drop from 0.4 
to 0.9. The increase of n(H®)/n(H) between adjacent cells is 
typically an order of magnitude when the larger (H®/H) is 
0.4 or more. 

We kept the number of cells within the model cube con¬ 
stant and varied the physical length of the cube proportional 

1 /3 

to ±49 . This makes clumped models with the same initial 
seed have the same ionization structure in terms of the size 
of the model cube. 

Fig. ID shows (a; upper left) the clumped gas densities 
in a typical slice through the plane containing the star; ( 6 ; 
upper right) the fraction of H® in the same slice, (c; lower 
left) the intensities of Hq, relative to its maximum intensity, 
as projected onto the sky (with no foreground extinction), 
and {d; lower right) the [NII]/Ha ratio projected onto the 
sky. The axes are in parsecs for a model with T* — 40 kK, 
I/49 = 10. We see that the density through the midplane is 
concentrated to the left side, but the projected images, (c) 
and (d), are more symmetrical because of ionization in the 
foreground and background planes. The [Nil]/Ha is large 
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Figure 1. Some properties of the models for L 49 = 10, T* = 40 
kK. The axes give distances in parsecs. Upper left; A slice through 
the midplane (containing the star), showing the densities in the 
models; upper right: the fraction of neutral H; lower left: a projec¬ 
tion of the Ha intensity upon the plane of the sky, normalized to 
the maximum intensity; lower right; the projected [NII]6583/Ha 
ratios. 


at the edge because the ionization favors N+ and the tem¬ 
perature is relatively high there. 

The large [NII]/Ha (lower right panel) at the edges does 
not imply that the HII region would appear limb-brightened 
in the sky. The lower left panel shows that the Ha intensity 
at the edges is only ~ 1 % as bright as the maximum, so the 
[NII] intensity is low and would blend in with the unrelated 
surrounding diffuse emission. It is only the ]NII]/Ha ratio 
that increases at the edges. 


3 EFFECTS OF CLUMPING ON IONIC 

ABUNDANCES 

Our results from a large number of models can be summa¬ 
rized simply: for a given exciting star, clumping changes the 
averaged fraction of ionization of all elements, but differently 
for relatively cool stars than for hot. The important quan¬ 
tity is 7 = (I/He/LH)/(He/H). If 7 > 1 a star is “hot”; < 1, 
“cool.” Clumping enhances the higher stages of ionization of 
an element (including He) for cool stars and decreases them 
for hot. The amounts of enhancement for specific ions are 
discussed below. 

We have Lue = f Vl/hvdv, integrated from the ion¬ 
ization potential of He'’, x(He°) = 24.6 eV, to x(He^) = 
54 eV. Lh is J Vf/hvdu, integrated from 13.6 eV to 54 eV. 


We count the radiation above 24.6 eV as both H- and He- 
ionizing because each recombination of He produces about 
one H-ionizing photon iOsterbrockll98flll . so absorptions by 
He hardly affect the number of H-ionizing photons. Inter¬ 
stellar dust does not affect the distribution of ionization sig- 
nihcantly if it contains only silicates, as seems likely (see 

3.1 Influences of parameters 

There are a large number of parameters in our models. Some 
set the basic physics and radiative transfer (elemental abun¬ 
dances, T*, L49, and /smooth), but others are arbitrary: the 
fractal dimension, D, and the initial seed that determines 
the density structure in clumped models. The dispersion of 
the results among the hve initial seeds that we arbitrarily 
chose for each clumped model is not large. Typical varia¬ 
tions of (H+/H) are 0.1%; of (He°/He), 0.2% if (He°/He) 
0.1 and <2% if (He'’/He)'~ 0.3. The mean temperatures 
weighted by nen(H''’) vary ~0.5%. Varying D from 2.6 to 
2.9 produced changes about twice the dispersion among the 
various seeds of each clumped model. The overall result is 
that the other parameters were much more influential than 
D and the initial seed. For the rest of the paper, we will 
present D = 2.6 results as representing clumped models, 
with results averaged among the five seeds. 

We begin by considering the main physical difference 
between clumped and uniform models, then dehne averaged 
abundances for comparing models, and then consider spe¬ 
cific elements and ions. 


3.2 Propagation of He Ionizing Radiation 

The most important consequence of a clumped density dis¬ 
tribution is the fate of He ionizing radiation, depending upon 
the hardness of the stellar radiation. It determines the rela¬ 
tive amounts of the nebular opacity supplied by H and He, 
the two big players for almost all ionizing photons. 

Hot stars produce bountiful He-ionizing photons, so H 
can compete with He for He-ionizing radiation. At hv = 25 
eV, just above x(He'’), H contributes ~40% of the opac¬ 
ity near the star with T* = 45 kK, and He practically all 
of the rest. In a uniform model, the fraction of H opacity 
gradually falls with distance from the star as He becomes 
more neutral. In the interclump regions of clumped models, 
the H opacity decreases less rapidly with distance because 
the density is lower. The abrupt increase of density when 
stellar photons encounter a clump makes the increase of He 
opacity over H occur within a small volume. The overall ef¬ 
fect is to make H compete for He-ionizing radiation more 
effectively in clumped models than in smooth, thereby de¬ 
creasing (He'''/He) and other high stages of ionization over 
uniform models. 

He supplies 95% of the opacity for 25 eV photons when 
T* = 37.5 kK. In this case, in clumped models the He is ion¬ 
ized in the interclump medium until the outflowing starlight 
encounters a clump, at which point both H and He tend to 
become neutral. In a uniform model, the He goes neutral 
after it absorbs essentially all of the stellar He-ionizing ra¬ 
diation. and there is a substantial region in which H is still 
H+ but He is neutral. The overall effect is that He and H 
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tend to be both ionized within the same volume in clumped 
models, so the (He+/He) within the ionized region is larger 
than in uniform models. The same effect applies to other 
high stages of ionization. This is the major physical cause of 
differences between clumpy and smooth models. 

We compared clumped models with those with uniform 
densities, but the root physical causes of the differences lie in 
the radial symmetry, not the uniformity, of the uniform mod¬ 
els. The general distinctions between clumped and uniform 
models applies to a ny smooth distr ibution of gas density. 

We considered iKuruczl (ll994^ stellar atmospheres of 
(37.5, 40, 42.5, and 45) kK, log g — i and used a power- 
law interpolation between them to estimate spectra for 38.1, 
38.7, and 41 kK as well. For these atmospheres, 7 = 1 oc¬ 
curs at about 39.1 kK. We will see qualitative differences in 
the ionization structure for other elements at this T*. For 
Kurucz atmospheres, 7 = (0.569, 1.35, 2.41) for T* = (37.5, 
40, 45) kK. We present plots of ionization fractions vs. T* 
because 7 , the real independent variable, is a less familiar 
quantity. 


3.3 Comparing Models with Observables 


Our aim is to compare the abundances within our clumped 
models with the abundances that would be inferred if the 
emissions from the models, measured without error or dust, 
were interpreted in the ways used on actual nebulae. The 
derived abundances differ from the actual ones in the model 
because of temperature and density variations within the 
model, even if the density is never enough to de-excite the 
emission lines by collisions. If there is poor spatial resolu¬ 
tion, there is averaging over the face of the nebula as well. 
We will compare three types of averaged abundances of a 
particular ion, X*, relative to H. The line emission (photons 
cm“® s“^ ster“^) is n{X^)ne ,Te), and jx is presumed 
to be known from atomic theory. 

Since the model provides 3-D ionization structure of the 
nebula, the true abundances are known: 

/n(He+)dK 

VH+Le /n(H+)dK ■ 

An observer would know the projected line intensities, 
Ix = f n{X'')ne jx{X'',Te) dz, where z is the distance along 
the line of sight, averaged over the observer’s spatial resolu¬ 
tion. An extreme assumption would be that there is no spa¬ 
tial resolution over the face of the nebula, so that we know 
only the total line fluxes, Fx, thereby obtaining a “global” 
averaged abundance: 


Fx 


H+ 


global 


J Ix dx dy , 

( 2 ) 

J n{X*)n,jx{X\T)dV . 

( 3 ) 

Fx{X^)/jx{X\Tx) 

( 4 ) 

F{H)/jH{TH) 

f n{XRw{X^)dV 

J n{H+)w{H)dV' 

( 5 ) 


where Tx is the measured averaged temperature relevant 
to the ion X*, depending upon whether it is a high or low 
stage. The weights w(X*) = He jx(X'',T), and similarly for 


the H lines. For high ions we use r([OIII]), obtained from 
[OIII]A4363/(5007+4859). For low ions we use r([NII]), 
from [NII]A5755/(6583-b6548). 

We will explore the effects of the ui(X*) on the derived 
abundances. Note that the effect is not limited to variations 
in temperature: it exists even in artificial isothermal models 
because of variations in n^. 

Another extreme assumption is that (X*/H) can be de¬ 
termined at each point (x, y), and the result, Ax{x, y), then 
averaged, weighting the bright pixels more than the faint. 
We define the abundances averaged over the face of the neb¬ 
ula by 

/ 4lx {x,y)lH{x,y) dx dy 
/ lH{x,y)dxdy 

The fluctuations of temperature within our models are 
modest, as uniform and smooth models have previously 
suggested. The variation of TOOI III) across the face of 
a nebula (e.g., lO’Dell ~ ZIIM is described by R = 
J(T(x,y) — ToRdx dy/lTg J dxdy], where dxdy is the ele- 
ment of area on the sky and To is the averaged temperature. 
lO’Dell et all (l2003ll found R in the range 0.005 - 0.016 in 
hve regions of the face of the Orion Nebula using HST im¬ 
ages of A4363 and 5007 with high spatial resolution. For T* 
= 40kK and L 49 = 10, our models collapsed along one axis 
yield R = 0.0012 for the uniform model and ( 1 . 0 — 1 . 6 ) x 10 “^ 
for the clumped models with various initial seeds. For T* = 
45 kK, R ~ 0.002. Thus, the Orion Nebula shows temper¬ 
ature fluctuations in excess of models, but not enough to 
explain the discrepancies between collisionally excited and 
recombination lines unless one assumes that in real nebulae 
more violent fluctuations are averaged out by the integration 
along the line of sight. 

We do not discuss fine structure lines ([O III] 52, 88 jam, 
plus several others) because their emissivities are not sen¬ 
sitive to variations in temperature or densities within the 
rather modest ranges that we consider. These lines provide 
excellent diagnostics of the true abundances of their respec¬ 
tive ions for both clumped and uniform models, but are sub¬ 
ject to collisional de-excitations in real objects. 

We now discuss individual ions and elemental abun¬ 
dances. 



3.4 He+/H+ 

Figure 121 displays (He'^/H+) plotted against T* for clumped 
models (upper sets of curves) and (He'''/H^ — 0 . 1 ) for uni¬ 
form density distributions (lower sets). The horizontal lines 
at (He^/H'*') = 0.1 and 0.09 show the cases of complete ion¬ 
ization of He. The solid lines are true abundances (eqn.0 for 
L 49 =0.1; the dashed lines are for L 49 = 100. The dot-dashed 
lines are the global abundances (that would be deduced by 
observers using the fluxes of lines) for L 49 = 0 . 1 , the dotted 
for L 49 = 100. The other values of L 49 have curves lying 
between the ones shown. The bends in the curves near T* 
= 39.5 kK occur because 7 ^ 1, so that the ratio of stellar 
ionizing photon luminosities for He and H, Lae/La equals 
the elemental abundance ratio (He/H). 

Fig. m shows that clumped models provide a larger 
He'''/H’'' at a given low T* than uniform models; ~0.06 at T* 
= 37.5 kK, as opposed to ~ 0.04 for uniform. For hot stars. 
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Figure 2. The ratio as determined for clumped (up¬ 

per set of curves) and uniform models, plotted against the excit¬ 
ing star temperature, T*. The uniform models show [(He'*"/H^)- 
0.01]. The horizontal lines show complete ionization of He for 
clumped and uniform models. Solid line are true abundances (see 
eqn. [3 for stars with L 49 = = 0 . 1 ; L* is the ioniz¬ 

ing photon luminosity. Dashes: true abundances for L 49 = 100. 
Dot-dashes: global clumped abundances that would be derived by 
observers (see eqn. ISJ for L 49 = 0 . 1 ; dotted, global abundances 
for L 49 = 100. The ratio of stellar ionizing photon luminosities for 
He and H, Lne/Tlt equals the elemental abundance ratio (He/H) 
at T* = 39.5 kK, producing the bends in the curves. 

the uniform models (the lower set of curves) have (He''"/H''') 
= (He/H) to high accuracy, but this is not so for clumped 
nebulae. They have ~2.5% of He'^ even for T* = 45 kK for 
the true abundances, and the derived global (He^/H'^’) (i.e., 
weighted by rie) is 98.1% of (He/H). The difference between 
the derived (He+/H+) and (He/H) is 3.8% for 42.5kK stars. 
The error increases for lower stellar temperatures and L 49 . 

The relatively large fraction of neutral He for clumpy 
nebulae surrounding hot stars has important consequences; 

(a) (He/H) is a major prediction of Big Bang Nu¬ 
cleosynthesis, and the comparison with real nebulae re- 
quires very high accura cy in the correction for He° (see 
lOlive fc Skillmanll 20041 for a review showing the accuracy 
required). The reason why clumped models exhibit more He° 
than uniform models is independent of metallicity. The ob¬ 
servations and interpretation of He I emission lines in metal- 
poor extragalactic H H regions strive for accuracies of ~1% 
in He^/H’^ (e-g-, IPeimbert et al.lr 2 nn(lll . There must be 
corrections for stellar absorption lines underlying the neb¬ 
ular emissions as well as density and temperature effects 
on the He I line strengths. Uniform or smooth models for 
nebulae ionized by collections of very hot stars predict that 
He''"/H''' = He/H. Clumped models suggest an increase of 
He/H of the order of ~ 3% to correct for He'^. Equation 11 of 
lOlive fc SkillmanI ll 20n4h suggests that uncertainties in the 
determination of the primordial He''"/H''" amount to ~10%, 
but the correction for He° is systematically upwards. 

(b) Helium is dominant in determining the opacity for 
energetic photons that produce ions of other elements, ex¬ 


ceeding hydrogen by a factor of a few near its ionization 
edge. For harder radiation {hv > x(0''") = 35 eV), its main 
competitor. O'*", typically provides <30% of the opacity, and 
the other ions contribute significantly less because of their 
relatively low elemental abundances. Thus, the increased 
He'’ decreases the abundances of high stages of ionization 
of heavy elements. 

3.5 Heavy Elements 

For heavy elements, the ionic abundances are strongly in¬ 
fluenced by whether He recombinations to the ground term, 
hiz ~ 25 — 26 eV, can ionize the stage in question (e.g., 
X(S+) = 23.3 eV), so that the He recombination radiation 
competes with stellar radiation. If the He radiation cannot 
produce ionization, for cool stars the lower stage (O^, Ne'*', 
N"*") is reduced in clumped models as compared to uniform. 
This behavior is similar to He” (figure because lower He” 
implies a smaller opacity for hard stellar photons. 

There are systematics in the temperatures measured in 
various ways. For all geometries, r([OHI]) is cooler than 
the rje-averaged temperature, (T) = J Tn^dVj f UedV, 
by ~300 K because is a powerful coolant. r([NH]) is 
warmer than (T) by ~400 K because stellar photons are 
always harder on average at the outer regions of the ion¬ 
ized gas where [N H] is produced, and is not as efficient 
a coolant as O^'^. Clumped models have warmer (T) than 
uniform; for T* = 37.5 kK, the difference is 50 - 100 K; for 
T* = 45 kK, it is ~ 180 K. r([NH]) for uniform models 
is about 50 K larger than for clumped. These temperature 
differences affect the derived abundances. 

3.5.1 High stages of Ionization: O’^'^ and Ne”'"'^ 

This section deals with the two ions with very similar be¬ 
havior as regards differences between uniform and clumped 
models: and Ne"*"'^. Their important characteristics are 

that they cannot be ionized by He recombination radiation 
and cannot be ionized to still higher stages by stellar radi¬ 
ation (x(0’'’’^) = 54.9 eV; x(Ne"''’^) = 63.4 eV). This group 
does not include (y = 34.8 eV) because it can be ion¬ 
ized to by stellar radiation. We will assume that there 
are only observations of optical lines and discuss nitro¬ 
gen in the next section, although there are fine-structure 
transitions of at 57 /rm and of at 69 pm. 

Figure El shows (O’^'^/O) plotted against T* for both 
true abundances and the global averages, for L49 = 10. 
Heavy lines are clumped models and thin are uniform. The 
true abundances are solid lines and the global abundances 
using r([OHI]) are dashed. 

Figure El showing (Ne''"'^/Ne), is analogous to Figure 
El The similarities in the ionizations are clear. Ne'*"'^ starts 
lower than at low T* but increases faster because it 
requires harder radiation: x(Ne''') = 40.96 eV; x(0’^) = 35.1 
eV. 

Figures El and El show that: 

1. For clumped models the accuracy of the global 
(0"’’’^/0) abundances derived from T([OHI]) is only ^70 - 
80% of the true. Even perfectly accurate [O III] and H/I line 
intensities do not lead to an accurate determination of 
abundances. The reason is that there are relatively strong 
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Figure 3. The fraction of oxygen in vs. T* for L 49 , the 

exciting stellar photon luminosity in units of 10 “^^ sec~^, = 10 . 
The heavy lines are clumped models; the thin lines, uniform. The 
solid lines are the true abundances; the dashed are derived from 
T([OIII]); the dotted are the spatially resolved abundances, av¬ 
eraged over the face of the nebula (see equation ( 3 ) 



Tyiooo K 

Figure 4. Plots of the (Ne'*"^/Ne) ratio, with the same nota¬ 
tion as used for figure 1^ Note the strong similarities, for reasons 
discussed in the text. 

temperature fluctuations within the zone in clumped 
models. Excess [O III] A4363 is produced by the warmer re¬ 
gions, making the mean r([OIII]) within the zone higher 
than the true mean temperature. This situation leads to 
an erroneously low abundance being determined for 
(Ne''"^/Ne) is somewhat more accurately determined than 

(o+yo). 

For uniform models, r([OIII]) global abundances are 
in excellent agreement with the true. Surprisingly, the ac¬ 
curacy of the r([OIII]) abundances is better for Ne”*"^ 



Tyiooo K 

Figure 5. Global (0”*"/0) vs. T* for L 49 = 10. The heavy lines 
are clumped models; the thin, uniform. The dashed lines are de¬ 
rived from r([OIII]); the dot-dashed from T([NII]). The dotted 
lines are averaged abundances. 


than for For both uniform and clumped models, 

(Ne+yNe)-( 0 +yO) until both are 

2. The averaged spatially resolved abundances (dotted 
lines) are almost the same as those obtained with no spatial 
resolution. 

3. The true (0’''^/0) and (Ne''"^/Ne) are larger for 
clumped models than for uniform in cool stars and smaller 
for hot stars, as was true for (He"*'/H"'’) (see Figure |^. 

The behavior of (0''"^/0) and (Ne"'"^/Ne) at other val¬ 
ues of L 49 is as expected. Decreasing L 49 decreases high 
stages of ionization: for L 49 = 0 . 1 , the set of curves has the 
same form as those already plotted but lie entirely under 
the envelope of the curves shown. 


3.5.2 Low Stages of Ionization: O’*' and N"*" 

These ions and their next higher stages have far larger abun¬ 
dances than the remaining stages of ionization (0°, N°, and 
over the bulk of the nebula. Both O"*' and N'*' are fur¬ 
ther ionized only by stellar radiation: = 35.12 eV 

and x(N’’") = 29.6 eV, well above the energy of almost all 
He recombinations. 

The (O’^/O) global ratios for L 49 = 10 are shown in 
Figure Heavy lines are clumped models, thin are uni¬ 
form. Solid lines are true abundances; dashed are global 
abundances derived from T([OHI]) and dot-dashed from 
T([NH]). The dotted lines are averaged abundances. Fig¬ 
ure (Ogives (N’''/N). Both show that the use of T([OHI]) 
for a low stage of ionization is problematic, especially for 
cool stars. The r([NH]) abundances for N”*" are ^40% too 
large but are, surprisingly, more accurate for ( 0 *^/ 0 ) than 
for (N+/N). 

Using T([OHI]) for and N’^ is especially problem¬ 
atic for clumped nebulae with hot exciting stars; the abun¬ 
dances will likely be overestimated by a factor of ~1.5. Us- 
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Figure 6 . Plots of the N”*" against T* for L 49 = 10 for clumped 
models (heavy lines) and uniform models (thin lines). Solid lines 
are true abundances. Dot-dashes: global abundances derived from 
T([NII]). Dashes: global abundances using T([Oni]). Dots: aver¬ 
age abundances (see equation |^. 


ing r([NII]) provides results within ~20% of the true for 
(0+/0) and ~30% for (N+/N). 

3.5.3 Neutrals: and N° 

In regions where these ions are abundant, their ratio to the 
singly-charged species is almost completely determined by 
charge exchange with H°. [01] A6300 emission is produced 
almost entirely in a physically thin zone at the outer edge of 
the H'*" zone. In this region there is appreciable H° (>0.5%) 
for charge exchange and a reasonably high Te, needed to 
produce the collisional excitation of [O I]. Newly ionized gas 
flows away from the neutral material, and gas dynamics, 
neglected in our models, plays an important role in the local 
temperature and ionization structure. Both the width of the 
zone affected by dynamics and the mean free path of ionizing 
photons are less than our cell size, so we resolve this region 
poorly. We do not present predictions for (O'^/O), which 
we find to be ~3% for L 49 = 0.1 and 1% for L 4 g = 100. 
These values are almost surely overestimates. Our (N°/N) 
is similarly unreliable. The cells in which 0.25^ (H'^/H)<0.9 
have a mass of ~8% of the highly ionized plasma in uniform 
models and ~30% in clumped. However, the temperature in 
the partially ionized gas is falling rapidly as H'*'/H decreases, 
so the emissions from the almost neutral gas are small. The 
region 0.9 < (H+)/H < 0.99 adds only ~1% to the mass 
and is quite cool (<6000 K), so it contributes a negligible 
amount of emission. 


3 . 5.4 Sulfur 

Sulfur has four relevant stages of ionization: S’*" (x = 
23.3 eV, so it is ionized by the most energetic He recom¬ 
bination radiation); (x =34.8 eV, similar to O'*'^); S’*"® 



Figure 7. Plots of the (S"*"/S) against T* for L 49 = 10 
for clumped models (heavy lines) and uniform models (thin 
lines). Solid lines are true abundances; dot-dashes, global abun¬ 
dances derived from T([NII]); dashes, global abundances using 
T([OIII]); dots, average abundances (see equation IHl. Compared 
with Figure!^ the same diagram for O"*", we see the influence of 
the ionization of S'*" by He recombination radiation that does not 
affect O’*". 



T 7 IOOO K 

Figure 8 . Plots of (S'^^/S) against T* for L 49 = 10 for clumped 
models (heavy lines) and uniform models (thin lines). Solid lines 
are true abundances; dashes, global abundances using T([OIII]); 
dots, average abundances (see equation O. 

(x =47.3 eV), and There is too little stellar radiation 
to make S"*"^ important in any of our models. 

Figures [7| and E show (S'’"/S) and (S"'"^/S) plotted 
against T*, both for L 49 = 10. The scheme of the lines is 
the same as in Figures These ions do not have the 

properties of the corresponding stages of oxygen (Figures 
and|3l. 
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Comparison of Figure |7| with Figure [^reveals that S"*" 
has a much smaller true fractional abundance than be¬ 
cause He recombinations ionize S’*'. (S^/S) varies far less 
with T* as well. A greater fraction of O is in O'*" than S in 
so the clumpy nature of clumped models makes more 
difference for S"^. The abundance of in clumped models 
is three or more times larger than for uniform. 

For clumped models, the global abundance of S'*" de¬ 
rived from T([OIII]) (the dashed lines in Fig. |7J is about 
twice the true abundances. The S'^ derived from r([NII]) is 
about 50% too large. 

FigurelSlshows the trends of (S^^/S). The bulk of sulfur 
is but both and decrease with increasing T*, 
showing the growing importance of 

The physically implausible assumption that all dielec- 
tronic recombination coefficients for S are zero leads to an 
appreciable shift of the abundances of S ions to the higher 
stages: from (S+: S+^: S+^: 8 +-*) = (0.172, 0.762, 0.066, 
5 X 10“®) for L 49 = 10, T* = 40 kK clumped models to 
(0.115: 0.702: 0.182: 0.003) with no dielectronic recombina¬ 
tion for S. Besides the obvious direct consequences for the S 
line strengths, the averaged electron temperature increases 
from 8751 K to 8812 K because is a major coolant. 
However, the effects of geometry, T*, and I /49 dominate the 
effects of the uncertainties in the dielectronic recombination 
coefficients of S. 

3.5.5 Total abundances: O, Ne, and S 

Assuming that we have measures of both r([OHI]) and 
r([NH]), how well can we estimate total elemental abun¬ 
dances from normally observed emission lines? Oxygen 
might provide the best case, since all of its ionization stages 
are observable. We include our estimates of O'^ in our discus¬ 
sion. In a real nebula it would be smaller than our prediction, 
while O’*" would be correspondingly larger. 

Figure]^ shows the estimated global abundances of O, 
Ne, and S for clumped models, relative to the true abun¬ 
dance. The estimates for oxygen indicate which temper¬ 
ature is used to estimate O'*". T([OHl]) is always used 
for The estimate for Ne is derived from scaling the 

Ne”*"^ abundance by the model predictions of (Ne/Ne”'"^). 
The line marked “S(no S’*"®)” simply ignores that ion. The 
S(ICF) is based on an em pirical Ionization Correc tion Factor 
(fCF) to account for S'*"® iKwitter fc Henrvl200il . ft relates 
S/(S+ -bS+^) to (0/0+): 

s = (s+ -b S+^) X ICF ; (7) 

log(/C'F’) = -0.17-b 0.18/3 - 0.11/3^-b 0.072/3® ; ( 8 ) 

/3 = log(0/0+) . (9) 

We see that the fCF provides an excellent estimate for the 
abundance of S. 

Figure irni is similar to Figure]^ except for uniform den¬ 
sities. For clumped models there are compensating errors 
when we use r([NH]) for 0 + and r([OHl]) for 0 +®: the 
O'*" is too large; the too small. For uniform models, 
both are only slightly too small. The net result is that the 
estimated abundance (the dot-dashed lines in figures ITHl and 
13 is found to be ~ 310 ppM for all T* and L 49 , both uni¬ 
form and clumped. Recall that in our models the true (0/H) 
= 330 ppM. 



T*/1000 K 

Figure 9. The total estimated global abundances of various el¬ 
ements. The elements and method of estimating its abundance 
are noted by the lines. For oxygen, “using T([OIII])” refers to 
the temperature used for O'*" and O^, and similarly for “using 
T([N11])”. T([OIIl]) is always used for 0+^. The S(ICF) is from 
an Ionization Correction Factor to account for S”*"® (see text). 



TyiOOO K 

Figure 10. The total estimated global abundances of various 
elements. The elements and method of estimating its abundance 
are noted by the lines. For oxygen, “using T([OIII])” refers to 
the temperature used for O'*' and O*', and similarly for “using 
T([N11])”. T([01I1]) is always used for 0+^. The S(ICF) is from 
an Ionization Correction Factor to account for S'*"® 


The importance of being able to use a suitable mean 
temperature [say, r([NH])] for low excitation ions is clear, 
especially for cool stars. If we are forced to use T([OHI]) 
for as well as O'*"'^, both estimates are almost always too 
large (especially for O'*"; see Fig. |3 dashes). They average 
about 400 ppM, and the worst (uniform, cool stars, L 49 = 
100) is a factor of 1.7 too large. The best case is for uniform 
geometry models excited by hot stars with L 49 = 100. Its 
predicted oxygen is almost correct, but the corresponding 
clumped models predict 20 % more than the correct value. 
Clumping introduces appreciable uncertainties. 

The plots of total abundances against T* for other val¬ 
ues of L 49 are quite similar to those shown in Figures |3 and 
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llOl for L 49 = 10. For L 49 = 0.1, the O using r([OIII]) for O"*" 
(the solid line) is relatively flat at ~ 120 % for both clumped 
and uniform models, an improvement over the situation at 
I /49 = 10. For 1/49 = 1, the curves are close to those in the 
figures. The total S without S”*"^ is also better for lower L 49 , 
since the levels of high stages of ionization increase with L 49 . 
For L 49 = 100, the total S without S'*"® varies from 110% of 
the true at T* = 37.5 kK to 83% at 45 kK for clumped and 
120 % to 66 % for uniform models. 


4 REMARKS AND SUMMARY 

The clumped ionization fractions differ significantly from 
uniform, and there are significant errors in estimating ionic 
abundances from emission line fluxes. (He"*'/H"*’) is signifi¬ 
cantly increased for cool stars by clumping, and (He'^/He) 
is appreciable (2 - 3%) for even very hot stars. Clumped 
models show less (O^^/O) and (Ne''"^/Ne) than uniform for 
hot stars and more for cool. For clumped, the abundances of 
both ions estimated by the standard methods using emission 
lines and the appropriate mean temperature, T([OIII]), is 
only ~70% of the true abundances. For uniform models, the 
derived ionic abundances agree well with true. Real nebulae 
probably show the same phenomenon, and no improvement 
in spatial resolution will resolve it. 

The strong effects of clumping on ionic abundances for 
a given stellar energy distribution (' H3.5II makes the determi¬ 
nation of the properties of the exciting star from the nebular 
ionic ratios even more problematic than it is with smooth 
models. Another uncertainty is the presence of carbonaceous 
dust, with a strong absorption at ~17 eV that makes 
the spectrum of the exciting star appear hotter because El- 
ionizing radiation is absorbed more readily than Ele-ionizing. 
Silicate dust makes little difference to El II regions because 
it has a low, almost wavelength-independent cross section. 

Real HII regions have irregular outer boundaries if they 
are bounded by ionization fronts. By means of models with 
clumped density distributions, we have investigated the ef¬ 
fects of the complex shapes of the outer boundary on the 
analysis of ionic and elemental abundances. We feel that 
our models contain the elements of a basic physical effect 
that pertains to real HII regions in all of their geometrical 
complexities: their boundaries are also complex in shape, 
more so than our m odels. _ 

We considered iKurucd (ll994^ stellar temperatures in 
the range 37.5 kK - 45 kK and ionizing photon luminosities 
of 10 "^® - s“^, all for an averaged density of 100 cm“®. 

Of course, other models of stellar atmospheres have been 
proposed. The most important parameter of any of them is 
7 , the ratio of He-ionizing photons to H-ionizing. We consid¬ 
ered 7 in the range 0.57 - 2.4. Our results apply rather well 
for densities n < 10 '* cm~® to objects with the same L^^jn 
as our models. With our limited spatial resolution we found 
no regions truly shadowed from the central star. Clumping 
makes significant differences on the true abundances of ions 
of N, O, Ne, and S, as well as those derived from emis¬ 
sion lines plus measured nebular temperatures, r([OIII]) 
and T([NII]). We considered two extreme limits of spatial 
resolution: either none, so that only the fluxes of the emis¬ 
sion lines were used, or averaging the abundances as de¬ 
termined pixel-by-pixel over the face of the nebula. There 


were no major differences introduced by the two methods of 
averaging. 

If r([NII]) is available for analyzing the low stages of 
ionization, the overall abundances of O, Ne, and S are rea¬ 
sonably accurate because of partially compensating errors. 
For cool stars (for which most Ne is Ne"*") in clumped mod¬ 
els, (Ne/H) is badly estimated from [Nelll] optical lines. 
All of these elements are well estimated for uniform models. 
However, using r([OIII]) for both low and high stages of 
ionization leads to overestimates of (O/H) of ~40 - 60% for 
both clumped and uniform models, for cool stars. 

The underlying physical cause of the differences be¬ 
tween clumped and uniform models is the flow of He re¬ 
combination radiation from zones of high ionization of H 
that can be in physical proximity to regions near the outer 
edge of the nebula, with relatively large amounts of H®. We 
discuss this effect in §3.1 and for each ion. 
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